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1 Introduction 



In recent years, pre-big bang (PBB) cosmology |] has established itself as a logically 
consistent alternative to standard slow-roll inflation ||. While the ultimate verdict 
on which, if either, of these two versions of cosmology will survive is left to future 
experiments, at the theoretical level each one has its own advantages and shortcom- 
ings. 

The trade-off follows from the fundamental physical difference between the two 
scenarios: in slow-roll inflation the evolution is from larger to smaller curvatures, 
while in dilaton-driven PBB inflation the opposite is true. Thus, while in the former 
case Planckian/string-scale curvatures are naturally reached at the very beginning, 
pushing the problem of initial conditions into the full quantum-gravity domain (see 
e.g. 0), in the latter that problem can be tackled within a perfectly known (albeit 
non-linear) set of equations (see e.g. 0). On the other hand, since large curvatures 
naturally occur at the end of dilaton-driven inflation[|, the exit problem becomes 
harder to analyse than in the conventional scenario. 

It is generally believed that the exit problem of string cosmology ||, i.e. achieving 
a smooth, non-singular transition from inflationary to FRW-type string cosmology so- 
lutions, should be solved by a combination of two kinds of effects, corresponding to 
two different classes of corrections to the low-energy, tree-level effective action. While 
higher-derivative corrections should be able to bound the growth of curvature || (or 
to allow a reinterpretation of the lowest-order curvature singularity), loop corrections 
0, accounting for the back-reaction on the geometry from particle production, should 
stop the growth of the coupling (of the dilaton) when a certain critical-energy con- 
dition |8j is reached. Recent proposals of cosmological entropy bounds || appear to 
support [jllj these ideas. 



Unfortunately, present techniques are still inadequate for giving a convincing an- 
swer to the question of whether any one of these desired effects does actually take 
place. The question of fully non-perturbative higher- derivative corrections amounts 



1 A phenomenologically interesting consequence of this difference is that, in slow-roll inflation, 
Planck/string-scale physics is "screened" from observation by the subsequent long inflationary era; 
in the PBB scenario, it is in principle accessible to observation. 



1 



to finding an exact (2D) conformal field theory. However, since supersymmetry is 
spontaneously broken in the case of cosmology, chances to get such a CFT in closed 
form are slim. 

Similarly, understanding the full effect of back-reaction would require knowledge 
of the full one-loop effective action up to four derivatives. In spite of brilliant work 



by Vilkovisky and collaborators [11] on this problem, these are still early times for 
applying the results so far obtained to our problem. Thus our aim here will be 
more modest: we will try to understand in which direction does the back-reaction 
from particle production work in PBB cosmology while it is still a small effect. No 
attempt will be therefore be made to achieve (a fortiori to describe) a complete exit. 



We thus took an easier road, recently opened by Iliopoulos et al. |T^] and by 



Abramo and Woodard [I3|, |]T4|. Their method was already checked, whenever pos- 
sible, against the more intuitive approach of Abramo, Brandenberger and Mukhanov 
15j . Apart from our use of a different gauge, known to be more appropriate in string 



cosmology, our work can be seen as a straightforward application of that method in 
the context of string cosmology. As a check of gauge invariance, we will also reproduce 
some results of [Q. 

Another (apparently strong) limitation of our work is that we consider loop cor- 
rections to solutions of the low-energy effective action, i.e. we do not include higher- 
derivative terms in the latter. Since the back-reaction effect turns out to be of the 
same order as a four-derivative term, it may look suspicious to neglect such terms 
in the tree-level action while keeping them in the correction. There is, however, a 
limit in which such an approximation is justified: it is a large- N limit, where N is 
the number of species of the produced particles. The one-loop back-reaction is of 
relative order g 2 Nl 2 R and thus dominates over a tree-level higher-derivative correc- 
tion 0(£ 2 S R) at sufficiently largeg g 2 N. Finally, it can also be shown that higher-loop 
effects on the metric, the dilaton and the antisymmetric tensor (axion) are paramet- 
rically smaller at large g 2 N and small g 2 , justifying the one-loop approximation. This 
observation makes it all the more urgent to devote further effort to a non-perturbative 
understanding of one-loop back-reaction in cosmology. 



2 N is indeed a large number for a typical string theory gauge group, say E s ® E s . 
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The paper is organized as follows: after recalling (Sections 2 and 3) the form 
of the lowest-order PBB background whose modification we wish to compute, and 
the general idea of the method, we will introduce (Section 4) our gauge choice and 
then go, in Section 5, to the heart of the calculation of the back-reaction, first from 
gravi-dilaton production, and then from gauge-bosons production. In Section 6, after 
expressing the final result in a convenient form, we will argue that the back-reaction is 
of the expected order of magnitude and that, at leading order at least, it modifies the 
original background in the right direction to help with the exit problem. Finally, an 
explicit check of the gauge independence of our results is presented in an Appendix. 

2 PBB Backgrounds and Fluctuations 

In units in which 167rG = h = c = 1 the normalized four-dimensional gravi-dilaton 
effective action in the Einstein frame reads: 



Even within this minimal set of moduli (corresponding to having frozen the axion 
as well as internal dimensions), the gravi-dilaton system offers several interesting 
cosmological solutions, including, of course, the simplest ones, which describe spatially 
flat, homogeneous and isotropic Universes. In the comoving (cosmic)-time frame the 
line element takes the form 



where t represents comoving time (we shall often use the suffix to denote tree-level 
quantities). In the discussion of quantum fluctuations it is often convenient to go over 
to the so-called conformal frame, in which the metric fl2.2|) becomes 



We recall that comoving and conformal times are related by dto = ttdr] . In the 
conformal frame a typical PBB-type solution appears as follows: 




(2.1) 
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where, by identifying the (arbitrary) constant rji with the beginning of dilaton-driven 
inflation (DDI), we have normalized Q to 1 at that moment. It is convenient to define, 
as usual, TC(rj) = Q'/Q = l/2rj, which is related to the physical Hubble parameter H 
by 

H« = ^ = % = H l ^ = ± 1 t <0, (2.5) 

do i V 61q 

where Hi = 1/27/ j refers to the initial value of the physical Hubble parameter. Note 
that, at initial time rj — rji, the physical Hubble parameter is given by Q' since Qi — 1. 

We remind the reader that the above background corresponds to DDI in the string 
frame, with a string-frame scale factor blowing up like (— t s )^ 1 /^ when the comoving 
string time t s approaches from negative values (t s — > 0~). In the Einstein frame, 
this is seen as a contraction (gravitational collapse, see @]). What remains frame- 
independent is the fact that curvature grows in time, rather than decreases as in 



power-law inflation. This is why the results of |L4| cannot be directly used in our 
context. 

Let us now define the fluctuations around the classical backgrounds by: 

g^ v = Q 2 (77^ + Vv) > (p = tp + 4> ■ (2.6) 

Our goal here is to compute the back-reaction on the above homogeneous backgrounds 
( |2.4| ) (i.e. ififu, and 4>) due to quantum fluctuations at leading order. Following flI3| , 
flU , we will thus parametrize the correction to the solution ( |2.4j ) by 

(4>) = D( V ) , (2.7) 

where t M = (—1, 0, 0, 0) is a time-like vector and fj^ u = rj^ + t^t u = diag (0, 1, 1, 1) . 
In the next section we shall recall the method of Refs. [12j , |13[ , |L4| for computing 
A, C and D. Here we just conclude by recalling that A, C and D have no physical 
meaning separately. For instance, if A and D are given in terms of C by 

'He 



D = f dr)'(nC)(r)') , (2.8) 

then all three can be gauged away by the redefinition of time, rj —>■ rj + J v QC. 
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3 The Procedure 



Since our background is of the power-law (though not of the power-inflation) type, 



we shall closely follow here the procedure described in [14] for computing A, C, D. 
We will then stress at which point the two procedures start to diverge. In |14| the 
authors first compute some amputated one-point functions denoted by a, 7 and S (see 
below), and then attach external leg propagators in order to reconstruct the desired 
functions A, C, D. Let us describe the procedure in a different, perhaps simpler, way. 

If, generically, we wish to compute a fluctuation over a generic background ^oi, 
then the relevant one-point function is given by the obvious functional integral 

f D^,e iS ^ oe+q,e ^(x) 

<*«<*» = jm^Jr ■ ^ 

If the action is expanded only up to the quadratic order in \E^, the one-point functions 
are trivially zero. If, however, the action is expanded up to cubic order, S^oe + 
#e] = Sl^o^ + S^l^oe^i^j + S^l^oil^i^j^k (the summation convention is adopted 
everywhere), and S^ k [^oi\ is treated perturbatively, then there is a non-trivial one- 
loop correction to the one-point functions (see Fig. 1). This is the familiar story 
of vacuum-shifting tadpoles in quantum field theory. Analytically the result can be 
expressed as follows: 

^ ~ J £^ /e iS[*M]-HS^[*M]*<** ' ^ ' ' 

The calculation of the diagrams in Fig. 1 can be performed in two steps. First, we 




Figure 1: The one-point functions from tadpoles 
compute the amputated one-point function for all possible three-point vertices, iS^ n . 
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Then the full one-point functions are obtained by attaching the appropriate external 



lines. Formula (|3.2|) gives 

(%(x)) = J dS(^(x)^(z)}(^ fe ( 2 )^ n (z)}[^f"] , (3.3) 

where (dropping indices hereafter), (ty(x)ty(z)) is the standard Feynman propagator 
for the quadratic action S^o]^ 2 = *A\I/, i.e. 

A(x)(%(x)y(z)) = i5\x - z) . (3.4) 

Finally, one can express the one-point function as: 

= -A- 1 (x)(^(x)^(x))[S 3 ] . (3.5) 

It turns out that the coincident limit propagator (\E'(x)\E'(x)) is a function of time, x°, 
only. Therefore, if one can invert the kinetic operator in the zero-momentum limit, 
then ( |3.5|) provides our final expression for the one-point function. For the purpose 
of the first step it is convenient to distinguish three kinds of amputated one-point 
functions, a, 7 and 5, defined respectively in Figs. 2, 3 and 4. 

As an illustration of the use of (33), we present here the calculation of the one- 
point function A in the covariant gauge introduced in Jl2). Note, from (|2~7) ), that A is 
nothing but the amplitude (ipn), which can be expanded from the general definition 
( [3. 3D and Fig. 1, in terms of the amputated one-point functions a, 7 and 5 as: 

a -> (M) ■ s + (M,*,)??' ■ 7 + (M^fT ■ « • (3-6) 



As in (|3.5| ), the external line propagators can be replaced by the inverse of the kinetic 
operators in the zero-momentum limit. In the covariant gauge of Iliopoulos et al. 
the Lagrangian can be diagonalized by transforming the variables into suitable linear 
combinations. The propagators of the original variables can be read off from those 
of the diagonalized variables by applying the inverse transformation. For a standard 
power- law-inflation background, Q ~ t s with s > 1, one easily finds 

s+1 s+1 

1 s 

(ihnl>w)1?r = —-iA A + —-iA c 
s+1 s+1 

(tMw>*T = f- 4 + -Ti) iA A + (3.7) 
V s + 1/ s + 1 
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where iAA and iAc are particular propagators of the diagonalizing variables. They 
can be replaced by the inverses of the corresponding kinetic operators in the zero- 
momentum limit, 1/Da and 1/Dc- With this replacement, when ( |3.7|) is plugged 
back into the expression for the one-point function A in ( |3.6|) , we immediately reach 
the desired expression for A: 



A 



1 

Da~ 



-4a + 



8 + 1 



(3a + j)-^-5 
s + 1 



+ 



D 



c 



S (3a + 7 ) + ^U 



s + 1 



s + 1 



(3.8) 



This equation is the same as was given in JL4J]. The above outline shows how simply 
the relevant formulae can be obtained once the general result ( p.5|) has been derived. 
In the Iliopoulos et al. gauge, which does need ghosts, there are three different types 
of loops, shown in Figs. 2, 3, and 4, where the dotted line describes the graviton, the 
dashed line the ghost, and the solid line the dilaton. 



••••• 



a = 



+ 



+ 




Figure 2: The a-type amputated one-point function 




Figure 3: The "f-type amputated one-point function 



Let us now briefly compare this approach to the one of Ref. JT^|, whose authors 
have used a more intuitive method by computing the energy-momentum tensor of the 
quantum fluctuations and their back-reaction over the original backgrounds directly. 
The quantum fluctuations are obtained by solving the linear order perturbation equa- 



tions and the constraints. In our approach, as in [13], one uses instead the constraints 
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- o - 



5 = j 



+ 



+ 




Figure 4: The 5 -type amputated one-point function 



and the gauge-fixing conditions to obtain the Feynman propagators entering the loops. 
In |15[ the back-reaction was computed through systematic corrections of Einstein's 
equations, the first non- vanishing correction being quadratic in the Planck length £p, 
the small expansion parameter characterizing quantum fluctuations to leading order. 
The counterpart to this in our approach is that the same parameter, £|>, controls the 
loop expansion in quantum gravity. Incidentally, this last observation neatly solves an 
apparent paradox: classical gravitational waves always carry positive energy while, as 
we shall see, the back-reaction from the quantum production of gravitational waves 
from the vacuum turns out to have, at least in some cases, negative energy. The 
point is that, in the case of the amplification of vacuum fluctuations, there is no pa- 
rameter allowing a separation of real-graviton production from other virtual one-loop 
effects whose sign is undetermined. This a-priori non-determination of the sign of the 
back-reaction is self-evident in the effective action approach ||11|| . 

We will follow the same procedure as outlined above in a different gauge where 
physical variables are more explicit. The usefulness of this gauge was first pointed 



out in |L6| and its significance has been explored in detail in fll7| , in the context of 
linear order perturbations to PBB backgrounds. In the next section we summarize 
results in the off-diagonal gauge relevant for the background ( [2.4|) . In the Appendix 
we will comment on the gauge dependence of the results by comparing them to those 
obtained in other gauges, e.g. in the covariant gauge of Jl2| . 
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4 Propagators in the off-diagonal gauge 



As explained in the previous section, in order to analyse the perturbations we need 
to expand Q2.1| ) to quadratic and cubic orders in the fluctuations. We find: 



S = J d A x^g g aP g pa g^ u \^ aP) ^vo,f3 - ^Vw.pVw + ^Vv,pVw 



(4.1) 



Here the 'pseudo-graviton' (Vv) indices are raised and lowered by 77^. To complete 
the expansion of the g^ v fields in terms of the background and the pseudo-graviton, 
we recall the expansions 



ft 4 



-g 



(4.2) 



where ip = 



Following [18[ we split the fluctuations into tensor and scalar parts by 




hij 



+ 



2$ -diB 
-d t B 2(*6 ij - didjE) 



(4.3) 



where represents the (transverse and traceless) tensor part with ha = and 
dihij = 0. The trace ip is given by ip = -2$ + 6^ - 2V 2 E, where V 2 = The 



action fl2.1|), when expanded up to quadratic order in the variables (|4.3j), shows that 
B and $ behave like Lagrange multipliers, giving rise to the following constraints (a 
prime refers to d/df]): 



$ : 4V 2 ^ - UH^' + 



4HV 2 (B - E') = . 



(4.4) 



Making use of the first constraint in (fO|), it is possible to diagonalize the quadratic 
part of the action in our background as 



S 2 = ± J d A x (v'v' - d iV d iV + ^v 2 ) + I(ci4 - d^iAQj + ^QjC 



(4.5) 
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where v and are the canonical variables, given in terms of the original fields as 



n 



1.1 



(4.6) 



Owing to two local gauge symmetries that preserve the form fl4.3|) , two out of the four 
functions can be fixed through two independent gauge choices. The off-diagonal gauge 
corresponds to choosing \& = E — 0, which simplifies the first canonical variable to 
v = Vtcf). The rest of the variables become related to this canonical variable in simple 
ways because of other relations following from the constraints (|4.4j) , $ = \/30/2 and 
V 2 5 = -V30/2. 

The process of obtaining the propagator for the canonical variables is simplified by 
dragging the scale factor, Q, from the canonical variables to the quadratic differential 
operator. In the off-diagonal gauge this amounts to 



S-2 



d A x 



^hijAhij 



A = Q 



-dl + V 2 



n. 



(4.7) 



The differential operator A, in the zero- momentum limit, is easily inverted: 



a 1 /(^) = - P d Vl n- 2 ( Vl ) r d V2 f( V2 ) , a = n 

Jlli Jr li 



-d 2 + — 



Q. 



(4.8) 



as has been advertised in |TJj]. This choice of the inverse operator also enforces the 
boundary conditions 



A^/fa) = (A - 7)%) = 0, 



(4.9) 



a property we require for the corrections A, C and D. 

From the quadratic action ( [4.7| ) it is easy to read out the Feynman propagators 
for the physical modes 



(0(x)0(y)> 
(hij(x)h k i(y)) 



iA 1 (x,y) 



2T, 



ijkl 



iA (x,y) 



(4.10) 



where is an appropriately normalized tensor enforcing the transverse and trace- 
less conditions on gravitons, i.e. diT ijk i = djT ijM = d k T ijk i = diT ijM = and 
T iik i = Tij kk = 0. For our purpose we are interested in computing the propagator 
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A 1 (x, y) in the coincident limit x — > y. To do this we go to the Fourier-transformed 
variables 

£^ 3 MvV k -\ h^x) = J ^^(^)e ik - x , (4.11) 

where both the mode functions 0k and h^j can be obtained by solving the differential 
equations A/k(r/) = 0. The solutions can be expressed in terms of Hankel functions 
of the second kind. Notice that, v = Q<f) being the canonical variable, one should 
normalize its modes in the far past, r\ — > — oo, as l/^/[k|. By the large-argument 
limit of the Hankel function, this completely fixes the modes: 

1 — T<D, 



h(v) = ^V^H^(\krj\) . (4.12) 
Consequently the propagator in the coincident-point limit takes the form 



lim^x)^)) = J dkk 2 H^(k\ V \)HP(k\ V \) 



lim (hi j(x)h k i(y)) = J dkk 2 (5 ik 5 j i + 5 a 5 jk - 5ij5 kl - S ik ^^ - bjr 2 ^- 

kjk k kiki kikj k k k\ kikjk k k\, 

~ dil n~r ~ d i k -gr + dkl ~pr + d ^-pr + — u — ) 

IX l\i IX, IXi f\i 

xH^(k\ V \)H^\k\ V \) . (4.13) 

Here we have replaced the transverse, traceless tensor by its Fourier transform. 
The integral in ( [4.13|) is in general divergent in the ultraviolet. It should be kept in 
mind, however, that the present approach should be complemented by a momentum 
cutoff as a result of the ultraviolet sickness of gravity. We set the cutoff, at any given 
time, at k max ~ which corresponds to a wavelength of the order of the Hubble 

radius. This is also the scale at which the mode functions switch from free plane waves 
to 'frozen out' modes. For modes larger than this wavelength the bosonic contribution 
dominates its fermionic counterpart (we assume some underlying supersymmetry) as 
a result of Bose-condensation. This is the effect we are after. For shorter-wavelength 
modes new degrees of freedom should appear, which would ultimately provide better 
ultraviolet behaviour, as in superstring theory. We also impose an infrared cutoff 
coming from a finite duration of dilaton-driven inflation, k m i n ~ \Hi\, corresponding 
to the maximal Hubble radius. In conclusion we restrict momenta in all loops to be in 
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the range: \Hi\ < k < Q\H\ = \Hi\/Q 2 . At this upper- momentum/short-wavelength 
cutoff, Q\H\ ~ 1/M, i-e. our cutoff corresponds to \krj\ ~ 1. Note that, in the 
off-diagonal gauge, as a result of the constraints, all fields can be expressed in terms 
of 0. Thus only one type of amputated loops, the 5-type, effectively arises on the 
right-hand side of (|3.5| ). 

5 Computation of the one-point functions 

In the first two subsections we shall compute the amputated one-point functions due 
to graviton, dilaton, and gauge-boson loops. In the last subsection the external zero- 
momentum propagators will be glued in to reconstruct the full one-point functions. 

5.1 Graviton and dilaton loops 

To get the three-point vertices we use our expansion of the action ( |4.1|) to cubic order 
in the fluctuations, S3. The results are summarized in tables 1 and 2, showing all 
the relevant three-point vertices. In table 3 we list all the coincident loops appearing 
in 5 and their respective contributions to the loop integrals (x = \kt]\). There is 
another way of doing the same computation, which involves, in general, four types of 
amputated loops in the off-diagonal gauge (Fig. 5). 



Figure 5: All possible types of amputated one-point functions in the off-diagonal gauge. 
Here the line represents all possible coincident propagators. 

The first type of loops, having hij in the external leg, turn out to be zero, since the 
loop (see both tables 1 and 2), in the coincident limit, always picks up a factor of Sij, 
which gives vanishing contributions because of the tracelessness of hy. The second 
type of loops, having in the external leg, also gives vanishing contributions, but 
for a different reason. Note that ^oi = diB can be integrated by parts to act on 
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the loop. The loops, however, being functions only of time in the coincident limit 
( |4.13|) give zero when hit by a space derivative. Thus only the third and fourth types 
of loops survive in the off-diagonal gauge, the 5-type and 7-type loops as defined in 
Figs. 3 and 4. 



1 


2 n l u 2 7 7 ifj , 1 p 1 ^L)i/j , 2 p 2 ^)ifj a3 /3 3 ^o)L 


2 


—H d% 3 n aia2 nP 2 P 1 'ib a (l)ib « (2)ib a (3^ 3 


Q 
O 


—ili0 2 rj tj yj ail3l {L)yj a2l 3 2 {Z)p a3l33 {6)t 


A 
^± 




5 


-fi 2 ^ 3 9 3 tl ^^/^3^ i/3i(1) ^ Q2/32( 2) V ; a3/33 (3) 


6 




7 


in 3 ^" ^^^77^^/3,(1)^/3.(2)^3/33(3) 


8 


|fi 2 Sr^ 3 r 7 ^^^Va 1 /3 1 (l)V'a 2 /3 2 (2)Va3/33(3) 


9 


Ifi 2 ^ 1 ^7 7 ^^«3 V , ai/3i ( 1 ) V , Q2/32 (2) V , a3/33 (3) 


10 


|fi 2 9 2 ai ^ 2 ^3^ 1 « 2 ^ i/3i(1)V;a2/32( 2)^ 3/33 (3) 


11 


ifi 2 ^^ 5 r/ Q4Q5 ^^^^^3/.3^ i/3i(1) ^ 2/32(2 )^ 3/33 (3) 


12 


ifi 2 9r^ 5 r/ Q4Q5 ^^^3/.3^ 2 ^^ i/3i(1) ^ 2/32(2 )^ 3/33 (3) 


13 


ifi 2 9 2 ai ^ 1 ^^r / ^3 V;ai/3i( i) V;a2/32( 2)^ 3/33 (3) 


14 


ifi 2 ^5 3 Q5 W5 ^^^-3^3^ V;Qi/3i(1) ^ 2/32( 2) V ; Q3/33 (3) 


15 


ifi 2 ar^ 5 r/ Q4Q5 ^^^-3^3^^ i/3i(1) ^ 2/32( 2)^ 3/33 (3) 


16 


ifi 2 ^ 1 ^ 1 ^-3 r/ /33/3 2V;ai/3i(1) ^ 2/32( 2)^ 3A( 3 ) 



Table 1: List of all Vv/SiV^^V >a 3 ftj vertices. While choosing an external leg one 
should permute 1, 2, 3 m eac/i vertex. 



1 


-f ^^ 6 r/^^0(l)^ 4a5 (2)^ 4/35 (3)t a6 


2 


f ^^ 4 ^^ 5 ^ 6/36 0(l)^ 5Q6 (2)^ 5/36 (3)t a4 


3 


v / 3^9r 5 r / a7a «0(l)^4« 5 (2)^ 7a8 (3)t a4 


4 


-2 v / 3^^ 5 r/ Q6 ^0(l)^ 4a6 (2)^ 5a7 (3)r 4 


5 
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Table 2: ^4// possible 4>i) ai i3 1 4>a2i32 an< ^ 4 >( t ) ' l l ) a 1 i3- l vertices. Again one should permute 
appropriate indices when considering each vertex. 
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(0,0') 




(0,0") 


-Jdxx A H^*{x)Ho W {x)/16H i 7T\ V \ 5 


(<9j0, 9j0) 


- J dxx 4 H^* (x)H^ 1) (x)/lQH i 7r\r ] \ 5 


((9j0, V>Oi) 


- v / 3/rfxx 3 ^ 1) *(x)^ (1) (a;)/32i7 i 7r|r ? | 4 


(<%0, Voi) 


v/3 / dxx A H ( v ] * (x)Hq W (x) /32H i TT \ rj | 5 




\/3 f dxx^H'^ix)^ 1) (x)/32H i n\r]\ 5 




-3 / dxx 3 H' {1) *(x)H'o W (x) /64^7r|r7| 4 


(djibni, dill)™) 


j u v / u v // * i /I 


\ IT VJ 7 IT VJ 1 


-3 / dxx A H' w * (x)H' (1) (x)/64H i 7T In 1 5 




3/do;a; 3 ^ 1) *(a;) J ffo (1) (^)/8i? i 7r|7/| 4 


(diij jk , dj^pik) 





{diijj jk , ditpjk) 


-3 / dxx 4 H^ 1} * (x)H^ (x) /8Hi7T\r]\ 5 


(0',0') 


-fdxx 4 H' {1) *(x)H' {1 \x)/16Hi<K\r]\ 5 




-3 / rfxx 4 ^ (1) *(x)/fo (1, (^)/8^7r|r/| 5 



Table 3: All possible coincident loop propagators appearing for the off-diagonal gauge. 
The prime in Hq 2 ^ refers to d/dx, otherwise to d/drj. 



Another important feature of the off-diagonal gauge is the absence of the Faddeev- 
Popov ghosts. This can be readily seen from the gauge transformations of the gauge- 
fixing functions \I> and E. Their gauge transformations have been given in 



V ^V + H£(x), E^E-£(x), (5.1) 

where £ and £ are the two gauge parameters associated with the coordinate trans- 
formations: rj —>■ rj — £ and x % — > x 1 — d^. Clearly, the Faddeev-Popov determinant 
is field-independent and hence the ghosts are decoupled, unlike what happens in the 
covariant gauge of (12| . 



As an example, we present a sample calculation, say the three-point vertex #5 
of table 2. For illustrative purposes, it is easier to make use of the constraints at 
the beginning (in contrast, in the computer code its more convenient to impose the 
constraints later). As we have argued, this does not make a difference, which is at 
the same time an important consistency check for our method. Making use of the 
constraints in the off-diagonal gauge presented in the previous section, the tr^^), 
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present in the vertex #5, can be replaced by the scalar fluctuation 0. As we said, in 
this case there is only one type of amputated loop S e s, that can be read off from the 
vertex #5: 

5 eS = ^[-^ V0'> + 4fli<#') + n 2 (<jxp)" + n 2 (d i( j>d i( j>)} . (5.2) 

Now, one should look at table 3, where various propagators have been listed in this 
gauge and pick up the leading-order terms from ( |5.2| ) as we approach the limit, r\ — > 
CP. Recall that Q 2 ~ \r]\, and hence all the terms in ( |5.2j ) are of the same order, 1/rf. 

5.2 Gauge-boson loops 

Let us now study the back- reaction produced by gauge fields to the background Q2.4Q . 
There are two distinct types of gauge fields appearing in the low-energy string effective 
action — the NS-type and the R-type. The NS-type couple to all the fundamental 
moduli while the R-type couple only to the metric. Let us first concentrate on the 
NS-type gauge bosons. In the canonical frame the action is: 

S em = \f d 4 xe^V=99" P 9 ua F^F p(7 , (5.3) 

which supplements the graviton-dilaton action ( |2.1| ) in the presence of the NS-type 
fields. Since the gauge-boson action ( |5.3j ) is already quadratic in the gauge fields, 
these decouple, to linear order, from other fluctuations. Also, interactions stemming 
from a possible non-Abelian gauge structure can be ignored. We shall work in the 
Coulomb gauge where the temporal component of the gauge field is zero and the 
spatial components satisfy <9jv4j = 0. Going over to a canonical gauge-field variable, 
the action can be written in the form 

S«° = ~ J d'xQAQ , A = -d 2 + V 2 - e^' 2 (e-^ 2 )" , (5.4) 

where the canonical field is given by Q — e~ v °l 2 Ai. Fourier expanding as before, the 
mode functions satisfy the homogeneous equation 

f d A* 

Ci = J 7^Ck,(^)e ik - x , AO = , k*& = , (5.5) 
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the solution of which can be written down [19| in terms of the Hankel functions of 
the second kind. Following the same normalization arguments as presented before, 
the Coulomb-gauge Feynman propagator for the canonical fields can be expressed in 
the form 

Mx)Ci(v)) = J (03i(^ - 1?) [V*vH?\kx )] [V^nH^(ky )} e^*-y) (5.6) 

where the index A = (V3 — l)/2. The cubic vertices have been listed in table 4 and 
the coincident loops are presented in table 5. 

For the R-type gauge bosons, the coupling to the dilaton being absent in the effec- 
tive action (|5.3|), the mode functions satisfy the free wave equations: Ck« + & 2 Cki = 0. 
Clearly, these modes are not parametrically amplified; and hence their contributions 
to the back-reaction are negligible. 

The relevant diagrams for gauge bosons are the same as in Fig. [5] except that here 
the solid line represents the gauge loop. Notice, from table 4, that all the three-point 
vertices are quadratic in the gauge fields and at the linear order there is no coupling 
between the gauge fields and the gravitational perturbations. Among the loops, see 
again Fig. ||, the ones that have hij and ipoi as external legs are zero for the same 
reason as presented before. The rest of the diagrams are of the 7 or 5 type. As an 
example, let us take the vertex #6 of table 4. It contributes to 5 c s as 

5 eS = Iexp(-po)[-(44) + (d k Aid k Ai)] . (5.7) 

The exponential factor is common to all the vertices and is ~ which cancels the 

inverse factor in the second coincident-limit propagator, as can be seen from table 
5. It requires some more effort to see that this is also happening for the first entry 
of table 5, but again, even for gauge bosons, the leading effect in the amputated 
one-point functions is ~ l/i] 4 . 
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1 




2 


4 eA Pl Y^) u 2 ^3 '/ ^ 1 /3 5 ^;^a 5 l,'->JV ; 04/34l 1 y | 


3 


| exp(-^ )df d 3 "V eAi A* 6 (2)A Pe (3)^ (1) 


4 


|exp(-^ )^ 5 9 3 a5 r ?a5ft r / ^^M a6 (2)^ 6 (3)^ 4/34 (l) 


5 


-exp(-^ )^ 5 ^ 4 ^A /35 (2)A Q5 (3)^ 4/34 (l) 


6 


iexp(-^ o )9r^ 6 r/ a5/36 ^««0(l)A /35 (2)A Q6 (3) 


7 


-±exp(-^ o )df^ 5 0(l)A ft (2)A Q5 (3) 



Table 4: possible 4>A ai Ap 1 and il) ai p x A a2 Ap 2 vertices. Again one should permute 
indices when considering each vertex. 





J dxx 2 [^{v/Vi)~^ /2 H { ^*(\kr]\)}' ■ 
[V=v(v/Vi)-^ /2 H?(\k V \)]'/8-K\ V \ 3 


(fiii -^-ii 1 A{ 2 ) 


J dxx A {t) /rii)-^H^* {x)H^ (x) /8vr |r/ 4 1 



Table 5: All coincident photon loops. All primes refer to d/drj. 



5.3 Adding external propagators 

Following ( |3.5|) we have first computed the amputated one-point functions in the 
off-diagonal gauge. As it turned out, after summing over all the vertices in tables 
1 and 2, the leading-order terms, as 77 — > 0, have the time dependence l/r/ 4 . We 
now attach the external line, i.e. the inverse of the zero-momentum kinetic operator, 
A . From (|4.8|) we note that our amputated quantities have to go through a double 
77-integral and multiplication by Q~ 2 . Thus the leading-order time dependence of 
the one-point functions is simply l/rf. Going through the procedure more carefully, 
we see that, apart for some numerical factors (which are different for different one- 
point functions), the leading-order contribution to all one-point functions is of order 
l/HiTj 3 . This result suggests that, in the context of PBB cosmology, the dominant 
modes are those near the upper momentum cut-off k ~ r]~ l . As we have already 
argued, contributions from modes beyond this cut-off, \krj\ ~ 1, have to be combined 
with those from other fields, e.g. from fermions. The resulting combined effect will be 
typical of a short-distance renormalization of the bare parameters of the Lagrangian, 
which should be mild in a supersymmetric theory. Apart from this, the physical 
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effect from short distances should be parametrically smaller than the effect we have 
computed. Considering Q2.4Q , Q2.5Q , and the formula Hi = 1/2^, the contribution we 
have found can be simply expressed, in our units, as: 

^ c < D ^w?~^= Hl - (5 ' 8) 

What about the ratio between the three one-point functions in the off-diagonal 
gauge? Our answer turns out to be simply: 

A=0, C = V3D, D = -A 1 5 cS . (5.9) 

We should point out that, while ( |5T9| ) appears to parametrize only certain forms of 
metric fluctuations, other contributions, such as (ipoi) or (ipij), i 7^ j, have been 
shown to be exactly zero. Also, the vanishing of a easily follows from the property 
ha = and from the coincident-point limit of the loop ( [4 . 1 3| ) . 

To reach ( |5.9| ) one could follow another path, which parallels more closely |13[ and 



nj] . One can take into account all types of loops -corresponding to all the variables 
appearing in the fluctuations ( |2.6| )- and use the constraints at the end. As a matter 
of fact, this also provides a good consistency check of the method. Here again, from 
very general considerations (cf. the three-point vertices in tables 1, 2 and 4), loops 
with ipoi and hij as external legs vanish, and only 5-type and the 7-type amputated 



functions survive. Thus the results ( |5.9|) follow in a different way: 

A = 0, C = V3D, D = -A" 1 (d + ^7) , (5.10) 



where (|5.10|) is identical to ( |5.9|) since the constraints provide the relation £ eff = 

We are now ready to discuss the physical implications of our results for the PBB 
scenario. 

6 Discussion 

The physical effects due to back-reaction are best analysed, in a PBB context, by 
going over to the string frame. In order to prepare for that, let us write down the 



18 



final outcome of our calculation in the off-diagonal gauge as 

(ds 2 ) = fi 2 [-(l - C)dr] 2 + dx- dx] = -dt 2 + a 2 (t)dx ■ dx 
(<p) = <f + D, (6.1) 

where C and D are as given explicitly in (|5.10|) . Since we limit ourselves to times for 
which C <C 1, the new comoving time t is related to the background comoving time 

by 



dt~n[l--C)dr} 



1 



-C)dt 



(6.2) 



Expressing the dilaton and the Hubble parameter H associated with the scale factor 
a in the new comoving time, and after using the relation ( |5.9| ) between C and D we 
find 



(<p) ~ --j= [ln(-t) - V3D 
H = A (lna) .i_[ 1 + V 3D 



(6.3) 



In order to go over to the string frame, we introduce the corresponding comoving time 
dt s = exp({(p)/2)dt and the scale factor a s = exp((^)/2)a. After a straightforward 
calculation the resulting string-frame Hubble parameter and dilaton can be written 

as 



(<P) 



1 



(y/3 + 1) [ln(-t s ) - (2^ - 3)D 



(6.4) 



It turns out that, in the Einstein frame, D(t) is a positive number times the square 
of the Hubble parameter. For the minimal gravi-dilaton system, we find: 



D(t) = 0.00A7H 2 (t) . 



(6.5) 



Although the coefficient of looks small, we have to recall that, in our units (167rG = 
h = c = 1), Planck's time is given by t 2 P = 1/16-7T. In other words, the (relative, 
dimensionless) correction to the tree-level background is: 



D = 0.2At 2 P H 2 , 



(6.6) 
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i.e. does not contain any particularly small number. 

When one goes to the string frame, in which the string length i s is a constant, the 



same result as (|6.6| ) for D follows, provided one interprets tp as the effective, time- 
dependent, Planck time tp = £ s e^^ 2 and H is expressed in terms of the string-frame 
Hubble parameter and of (now the dot refers to derivative w.r.t. cosmic string 
time) by the standard relation H Q = H s — 0/2. 

When the NS-gauge-boson loops are included the numbers in ( |6.5| ) change but the 



conclusions do not. The contribution of a single gauge boson to D turns out to be 

D = 0.00016# 2 = 0.008 t 2 P Hl , (6.7) 

i.e. of the same sign but about 30 times smaller than that of the gravi-dilaton system 
(|6.5D. The result ( p-7\) should be multiplied by the number of NS gauge bosons, which, 
in typical superstring theories, could be in the range of hundreds. 

In both cases the one-loop correction adds to the rate of inflation and to the rate 
of growth of the dilaton. At first sight, this seems to go against an eventual exit from 
inflation. However, what is crucial for the exit is not so much the absolute rates of 
growth but rather the relative rate of dilaton and scale-factor growth. 

To be more precise, let us compute directly the effect of back-reaction on the 
string-frame quantity Tp, where Tp = tp — 3 In a s is the duality- invariant shifted dilaton 
and, as before, the dot refers to the time derivative in the string frame. From Eq. 
^4]) we easily compute: 

9-3^3 



1 

<p= — 

f s 



-D 



(6.8) 



The positive sign of D is forcing Tp to decrease (relative to the tree-level value) as the 
result of a competition between a growing dilaton and a growing Hubble parameter in 
the string frame. It is thus obvious that the ratio Tp/ H s will change from its tree-level 
value of \/3 towards smaller values. As shown in Fig. ^|, this is precisely as required 
by a graceful exit ||, i.e. the phase-space trajectory in the Tp-H plane should turn 
anticlockwise. Incidentally, this is also what the Hubble entropy bound of || requires 



1C 



Another way the effect (|6.8|) can be understood is by computing the effective 
energy and pressure densities produced by the quantum fluctuations and by check- 
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Figure 6: The phase diagram of PBB cosmology. The curve shows the direction in 
which the back-reaction we computed is modifying the lowest order solution (straight 
line) helping an eventual flow from DDI to a FRW-phase. 

ing whether the condition Sp s + Sp s < is obeyed, where Sp s and 5p s are the en- 
ergy and pressure densities in the quantum fluctuations. Let us study this first 
in the Einstein frame. For an additional ideal fluid of energy-momentum tensor 
T\f = diag(— Sp, Sp, Sp, Sp) added to the dilaton the energy density and pressure den- 
sity appear in the Einstein-Friedmann equations as 

Sp = 6H 2 - ^y? 2 , -bp = AH + 6H 2 + ^ 2 . (6.9) 

When both Sp and Sp vanish, we recover, of course, the tree-level solutions. Inserting 
instead the corrected metric and dilaton backgrounds, we easily obtain: 

5p = ~tS ) 6p = ' (6 ' 10) 

Since D is positive, see (|6.5| ) and ( |6.6|) , the energy density is negative in the Einstein 
frame. The value Sp = is special to the background Q2.4|) . We have checked that, for 
other power-law backgrounds (e.g. for radiation-like solutions) Sp does not vanish. 
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To translate these results to the string frame, we note that the corresponding 
energy and pressure densities are given by 



e*8p a = e^Sp = -3H 2 S + Tp 

e^6p s = e^dp = 2H S - 3H 2 S - 2H s xp - ip 2 + 2xp , (6.11) 

so that, again, 5p s < and 5p s = 0. Hence, the condition 5p s + 5p s < is clearly 
met, see Fig. |6|. 

We conclude that, from all points of view, the computed back- reaction goes in the 
right direction to favour a change to the FRW branch and to avoid violation of the 
entropy bounds. 

Appendix 

Here we compare our results with similar calculations done in other gauges, such as 
the Iliopoulos et al. gauge, which has been extensively used in the recent literature to 
compute back- reaction to chaotic ]13[ and power-law inflation |14| . Specifically, in the 



latter case, the authors found no strong effect (i.e. growing with time) at one loop. 
It might not be so clear to the reader why the one-loop effect is so dominant in PBB 
inflation, which is also essentially power law. To clarify the situation it is important 
to recall the basic distinction between PBB and standard power-law inflation. In 
the PBB case the Hubble parameter grows with time, as t — > 0~, and as a result 
higher frequency modes leave the horizon at higher values of the Hubble parameter, 
see Fig. [7]. The loop integrals are thus dominated by the latest value of the Hubble 
parameter, in contrast with the case of standard power law, where they are dominated 
by its initial value. It is also important to keep in mind the two different time limits: 
in the case of standard power law, the limit is t — > oo, whereas, in PBB it is, t — ► 0~. 
This makes a crucial difference when deciding which are the dominant or sub dominant 
effects. 

As a check of gauge invariance, we feel it is worthwhile to redo the standard 
power-law calculation in the off-diagonal gauge. It will provide a direct check of the 



results obtained in the gauge used in |14fl , as well as help to better understand the 



differences between the two cases in the same gauge. Let us recall some of the basic 
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Figure 7: Horizon crossings of various scales in PBB cosmology. 



formulas for standard power-law inflation. The scale factor, Hubble radius, etc., are 
given by (s > 1): 

lna = slnt = ^—hiT], H = H { a 1/s = , Hi = — -— (6.12) 

S - 1 t S - 1 7fc 

and the dilaton is given by (p = 2^/~s~ lnt . For an arbitrary value of s, ( |6.12|) describes 
a background solution only in the presence of a potential 

Vfr) = 2Hf (3 - i) exp [-<p/Vs\ . (6.13) 

Note that for the special value s = 1/3, corresponding to the PBB case, the potential 
vanishes. To proceed as before, in the off-diagonal gauge, we have to first find the 
quadratic constraints through which the metric perturbations become related to the 
canonical variable 0. One finds 

•=-2^*' V2B =2V- (6 ' 14) 
The rest of the analysis is almost the same as before, except one should take into 
account the extra three-point vertices coming from the potential (|6.13|) . Using the 
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constraints and the gauge choice, the one-point functions turn out to have the follow- 
ing forms 



where 6 and 7 represent the same amputated diagrams as before. After summing over 
the vertices, it turns out that the leading effect in the one-point functions is of the 
formf] D ~ lnt - If we proceed, as in the case of PBB, to study the contribution of C 
and D to ip and if one -ioop, we find that the leading logarithms disappear, irrespective 
of the numerical factor in front of Into- The reason for the cancellation is twofold: 
firstly, the relation (|6.15| ) between C and D insures that the changes in ip and H 
are proportional. Secondly, it is easy to show that for C ~ Into the change in H is 
subleading since Ct/C ~ (1/lnt) — > 0. 

Finally let us briefly point out another technical difference between the PBB and 
the standard power-law cases. In the PBB case the constraints relating $ and B to 
have signs opposite to those in ( |6.14j ) implying that C and D have the same sign, see 
( |5.9|) . This is why their effects on ip and H one _\ oop add up, in contrast to the case of 
standard power law, where they compete with each other. Furthermore, the leading 
behaviour is no longer logarithmic, explaining why a large effect survives only in the 
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(6.15) 



PBB case. 



3 For instance, for the particular case s — 2, we found D = 
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